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The effect of electric and color charge neutrality on 2-flavor color super- 
conductivity has been investigated. It has been found that the effect of the 



charge neutrality on a 2-flavor quark system is very different from that on 



a 3-flavor system. The BCS diquark pair in 2-flavor color superconducting 
phase has been largely reduced by the large Fermi momentum of electron, the 
diquark condensate first increases with the increase of quark's chemical poten- 
tial fi, then decreases rapidly and finally disappears at about /x = 535MeV, 
at which the thermodynamic potential equals to that in the neutral normal 
2-flavor quark matter. 
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I. INTRODUCTION 



Because of the existence of an attractive interaction in the anti-triplet quark-quark chan- 
nel in QCD, the cold and dense quark matter has been believed to favor the formation of 
diquark condensate and being in the superconducting phase [!]]- M. 

The QCD phase structure at high baryon density is determined by how many kinds of 
quarks can participate in the pairing. If we have an ideal system with two massless quarks 
u and d, and the three colored u quarks and three colored d quarks have the same Fermi 
surface, the system will be in 2 flavor color superconducting phase (2SC); and if we have an 
ideal system with three massless quarks u, d and s, and all the nine quarks have the same 
Fermi surface, then all the quarks will participate in pairing and form a color-flavor-locking 
phase (CFL) §. 

In the real world, the strange quark mass m s should be considered, which will reduce 
the strange quark's Fermi surface. If m s is large enough, there will be no pairing between 
us and ds, and a 2SC+S phase will be favored 0. 

Also, in reality, the system should be electric and color neutral, the existence of eletrons 
in the system will shift the Fermi surface of the two pairing quarks, and the system can be 
in BCS phase or crystalline phase J7| or normal quark matter, depending on how large the 
difference in the chemical potentials of the two pairing quarks is. 

In a recent paper ||, it has been argued that the 2SC+S phase would not appear if the 
electric and color charge neutrality condition is added, which is also agreed with the results 
in |§ by using SU(3) NJL model and taking into account the self-consistent effective quark 
mass rrif(fi), where / refers to u, d, s. It is found in || that the charge neutrality has a large 
effect on the s quark mass, i.e., m s (fi) begins to decrease at a smaller fi ~ AOOMeV than the 



case if no charge neutrality is considered like in |T0[ [|TTJ, where m s (n) begins to decrease at 
about /I ~ 550MeV . 

In this paper, complementary to ||, we investigate the effect of charge neutrality on 
2-flavor color superconductivity assuming that there is no strange quark involved in the 



chemical potential regime \i < 550MeV. 



This paper is organized as following: in Sec. II, we extend our method in |12j to derive 
the thermodynamic potential when charge neutrality is considered; the gap equations and 
charge neutrality conditions will be derived in Sec. Ill, and in Sec. IV, we will give our 
numerical results; the conclusions and discussions will be given at the end. 



II. THERMODYNAMIC POTENTIAL 



A. The Lagrangian 



Assuming that the strange quark does not appear in the system, we use the SU(2) 
Nambu-Jona-Lasinio model, and only consider scalar, pseudoscalar mesons and scalar di- 
quark. The Lagrangian density has the form as 

C = q(i^d^ - m )q + G s [{qqf + {qi^rq) 2 } + G D [{if ee b l5 q){iqee h lb q C )l (1) 

where q c = Cq T , q° = q T C are charge-conjugate spinors, C = i r y 2r )° is the charge con- 
jugation matrix (the superscript T denotes the transposition operation), the quark field 
q = q ia with i = 1, 2 and a = 1,2,3 is a flavor doublet and color triplet, as well as a 
four- component Dirac spinor, r = (i -1 ,-/ -2 ,-/- 3 ) are Pauli matrices in the flavor space, and 
(e) lk = e lh , (e b ) a/3 = e af3b are antisymmetric tensors in the flavor and color spaces respec- 
tively. Gs and Gd are independent effective coupling constants in the scalar quark- ant iquark 
and scalar diquark channel. 

After bosonization, one can obtain the linearized version of the model (|I|) 

£ = q{i^d^ - m)q - \/\*\if ee b ^q) - \& b {iqe€ b 1 sq C '" 
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where we have assumed that there will be no pion condensate and introduced the constituent 
quark mass 



m = m + o. (3) 



From general considerations, there should be eight scalar diquark condensates [|13j |n 
In the case of the NJL type model, the diquark condensates related to momentum vanish, 
and there is only one + diquark gap with Dirac structure T = 75 for massless quark, and 
another + diquark condensate with Dirac structure T = 7075 at nonzero quark mass. In 
this paper, we assume that the contribution of the diquark condensate with T = 7075 is 
small, and only consider the diquark condensate with T = 75. The diquark condensate with 
Dirac structure 7570 has been recently discussed in |L5] . 

We can choose the diquark condenses in the third color direction, i.e., only the first two 
colors participate in the condensate, while the third one does not. 

The model is non-renormalizable, and a momentum cut-off A should be introduced. The 
parameters G$ and A in the chiral limit m = are fixed as 

G s = 5.0163GeV~ 2 , A = 0.6533GeV. (4) 

The corresponding effective mass m = 0.314GeV, and we will choose Gd = 3/4Gs in our 
numerical calculations. 

B. Partition function and thermodynamic potential 

The partition function of the grand canonical ensemble can be evaluated by using stan- 
dard method Jl6| fl7|l , 

Z = N'J [dq\ [dq]exp{J^ dr J d 3 x (£ + nqioq)}, (5) 



where j3 = 1/T is the inverse of temperature T, and fi is the chemical potential. When 
electric and color charge neutrality is considered, the chemical potential fi is a diagonal 6x6 
matrix in flavor and color space, and can be expressed as 

fi = diag(n u // 2 , /i3, A*4, /x 5 , He), (6) 



A 



where fix, fi2, 1^3 are for the three colored u quarks, and ^,^5, A*6 are f° r the three colored d 
quarks. 

Like in m ||, the chemical potential for each color and flavor quark is specified by its 
electric and color charges fa = /jl — Q e fi>e + T 3 fi 3c + T 8 /z 8c , where Q e , T 3 and T 8 are generators 
of [/(1)q, £7(1)3 and U(l)$. Because the diquark condenses in the third color direction, and 
the first two colored quarks are degenerate, we can assume fi 3c = 0. For the same flavor, the 
difference of chemical potentials between the first two colored quarks and the third colored 
quark is induced by /x 8c , and for the same color, the difference of chemical potentials between 
u and d is induced by fi e . 

The explicit expressions for each color and flavor quark's chemical potential are: 

2 1 

Hx = 11-2 = n - -/M e + -/i 8c , 

1 1 

/i4 — A*5 — ^ + g^e + ~^8c, 

2 2 

^3 = H ~ T^e ~ -//8c, 

1 2 ,»\ 

06 = V + g0e - g/^8c- (7) 

For the convenience of calculations, we define the mean chemical potential p, for the 
pairing quarks qx,Q5, and 52, <Z4 

fiX+ fi 5 A*2 + 04 1 1 , Q x 

A* = 2 = 2 = ^ ~ 6^ e S^ 8 "' 

and the difference of the chemical potential Sfi 

r 05 _ 01 04 ~ 02 /o / n x 

<$0 = 7, = 2 = ^ 

Because the third colored u and d, i.e., the 3rd and 6th quarks do not participate in the 
diquark condensate, the partition function can be written as a product of three parts, 

2 = Z const Z^Zx^^i- (10) 

The constant part is 



Zconst = N'ex P {- f dr I d 3 x + ^}. (11) 
Jo J 4G5 4Gd 

Z 3e part is for the impairing quarks 53(1*3) and (ft (^3), and ^15,24 part is for the quarks 
participating in pairing, q\ {ui) paired with q 5 (<i 2 ), and ^2(^2) paired with g 4 (di). In the 
following two subsections, we will derive the contributions of Z 36 and Z 15 24 . 



1. Calculation of Z; 



36 



Introducing the 8-spinors for q 3 and q§, 



^36 = (93, ?e; ?f , 



(12) 



we can express Z 36 as 



'36 



Det 1/2 (^[Go r 



36 J 



(13) 



where the determinantal operation Det is to be carried out over the Dirac, color, flavor and 
the momentum-frequency space, and [Gq ^36 has the form of 



000^ 



[G 



t\* x 



[GolJ 1 
[Go], 1 



with 



(14) 



(15) 



Here we have used $ = p^7 M and fa = /^7 . 



For the two quarks not participating in the diquark condensate, from Eq. (|T3|) , we can 
have 



lnZ 36 = -lnDetO^Go 1 ^) 

= iln[Det(/3[G V)Det(^ (16) 

We first work out 

[Det(/?[G + ] 3 - 1 )Det(/3[G -] 3 ')}= P'[pI - EfM - E%\ 
[DetC9[G3-]^)Det(^[Go]5- 1 )]= ?M ~ EfM - Ef], (17) 

with Ef = E ± and E^ = E ± \x§ where E = \/p 2 + m' 2 . Considering the determinant 
in the flavor, color, spin spaces and momentum-frequency space, we get the expression 

n f 

+H(3 2 [pl - Ef]) + Mf[pl - Ef])}. (18) 



2. Calculation of Z\§ 24 



The calculation of ^15,24 here is much more complicated than that when the two pairing 



quarks have the same Fermi surface [12 



Also, we introduce the Nambu-Gokov formalism for q 1 , q 2 , q± and q 5 , i.e., 



* = Q2, 94, 95; q?, q%, q% , q°)- 



(19) 



The ^15,24 can have the simple form as 

^15,24 = /[d*]exp{±X; 



-1 



15, 24 J, 



(20) 



where 



P -i 

^15,24 



^ [^o] 15*24 A ^ 



(21) 



with 



[^0 ] 15,24 



([Gt]? X 

[Gt^ 1 

[G^I 1 

[G^ 1 



and the matrix form for A 1 * 1 is 



A" = -zA 75 



From Eq. fl20|) , we have 



^0 1' 



0-10 
0-100 
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, A + = 7 °(A-)V. 



lnZ 15i24 = -lnDet(/?G^ 24 )- 



For a 2 x 2 matrix with elements A, B, C and D, we have the identity 

Det = Det(-OB + CAC~ l D) = Det(-5C + BDB A). 

\cd) 

Replacing A, B, C and D with the corresponding elements of G^ 2 ^ we have 



Det(G^ 24 ) = DetA = Det(— A + A~ + A^^A-]" 1 ^]^) 
= DetD 2 = Det(-A-A + + A-[G^ 5 ) 24 [A+]- l [G+}^ 24 ). 



Using the massive energy projectors A± in [T2|] for each flavor and color quark, we 
out the diagocal matrix D\ and D 2 as 



(Di)u 


= Pl)22 


= [(P0 


+ <5p) 2 




+ [(po + 5/i) 2 


- [^] 2 ]A + 


(A)33 


= Pl)44 


= [(Pb 


-5/i) 2 


- [^aHA- 


+ [(po - <5p) 2 


- [^] 2 ]A + 


P 2 )ll 


= (£2)22 


= [(Pb 


-5/i) 2 


- [^] 2 ]A- 


+ [(po - <5p) 2 


- [^a] 2 ]A + 


(^2)33 = 


= = 


= [(Po H 




- [^] 2 ]A- 4 


- [(po + 5p) 2 - 


- [^a] 2 ]A + , 



where = \J{E± p) 2 + A 2 . 

With the above equations, Eq. (E3p can be expressed as 

lnZ 15i24 = 2 £ E{ln[/? 2 (^ - + <V)) 2 ] + M/^o - - S^f] 

n p 

+ \n[(3 2 (p 2 - {Et + 5/i)) 2 ] + \n[f3\pl - (£+ - 5/i)) 2 ]}. (28) 



C. The thermodynamic potential 

Using the helpful relation 



n 

= ^ p + 2Tln(l + e-^)], (29) 



we can evaluate the thermodynamic potential of the quark system 

,2 A2 



+Tln(l + exp{-(3E^)) + Tln(l + exp(-/3E£)) + Tln(l + exp{-(3E 6 )) 

+2EX + 2E A + 2Tln(l + exp(-pE+ + )) + 2Tln(l + exp(-(3E+^)) 

+2Tln(l + exp(-/3£ A+ )) + 2Tln(l + exp(-/3E^))] (30) 

with j£± ± = E^± 5fi. 

For the total thermodynamic potential, we should include the contribution from the 
electron gas, fl e . Assuming the electron's mass is zero, we have 

«. = -4- 2 . Pi) 

The total thermodynamic potential of the system is 

Q = Q q + Q e . (32) 



III. GAP EQUATIONS AND CHARGE NEUTRALITY CONDITION 



From the thermodynamic potential Eg. (p2|) , we can derive the gap equations of the order 
parameters m and A for the chiral and color superconducting phase transitions. 

A. Gap equation for quark mass 

The gap equation for quark mass can be derived by using 

dm ^ 

The explicit expression for the above equation is 

E^ 



-/(£++)-/(£+_)) 
+ (2 - f(E+) - f(Ez) - f(E+) - f(E 6 ))} = 0, (34) 

with the Fermi distribution function f(x) = 1/ \exp{/3x} + 1). m = corresponds to the 
chiral symmetric phase, m^O corresponds to the chiral symmetry breaking phase. 



B. Gap equation for diquark condensate 

Using 

we can derive the gap equation for diquark condensate 

A[l - 4G D / - REl + ) - f{Bl-)) 

+2^(1 - f(E£ + ) - /(£+_))] = 0. (36) 



i n 



C. Color charge neutrality 



The color charge neutrality condition is to choose /isc such that the system has zero net 
charge T 8 , i.e., 

T S = = 0. (37) 
Evaluating the above equation, we have the color charge neutrality condition as 

- f(EX + ) - /(£+_)) 
+(f(E+) - f(E-)) + (/(£+) - f(E-))} = 0. (38) 



D. Electric charge neutrality 

Similarly, the electric charge neutrality condition is to choose \i e such that the system 
has zero net eletric charge Q e , i.e., 

Qe = ^ = 0. (39) 



From the above equation, we obtain 
d 3 p ,£ 



/^[2^(l-/fe)-/fc)) 
-2^(l-M + )-M-)) 



El 



+6[f(E- + ) + f(E+ + ) - /(£+_) - f(E~_) 



+4(f(Et) - f(E-)) - 2{f{E+) - j{E-))\ + £| = 0. (40) 



IV. NUMERICAL RESULTS 



In the numerical results, we investigate the following four different cases: 1), no color 
superconducting phase, no charge neutrality condition; 2), no color superconducting phase, 
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electirc charge neutrality condition added, the electron's chemical potential in this phase 
will be characterized by 3), color superconducting phase exists, no charge neutrality 
condition, the diquark condensate in this phase will be characterized by A ; 4), color super- 
conducting phase exists, electric and color charge neutrality condition added, the electron's 
chemical potential and diquark condensate in this phase will be characterized by \i e and A 
respectively. 

Fig. m is for the chiral phase transition in the case of no electric neutrality considered, 
which is familiar to all of us, the quark posses a large constituent mass in the low baryon 
density regime due to chiral symmetry breaking, and restores chiral symmetry at high baryon 
density. 

Fig. |2] is the chiral phase diagram when the electric charge neutrality is considered, the 
solid squares and the solid circles are for quark mass m and electron's chemical potential yiP e 
respectively. It is found that, in the chiral symmetric phase, there is a large Fermi surface 
for electrons, and the electron's chemical potential /x° increases linearly with the increase of 
quark's chemical potential \i. 

Fig. H| is the phase diagram for color superconductivity without charge neutrality con- 
straints, the quark mass m (squares) and diquark condensate A (circles) are plotted as 
functions of quark's chemical potential /i. It is found that in the color superconducting 
phase, the magnitude of the diquark condensate A is about lOOMeV by using the param- 
eter Gd = 3/ AGs- A first increases with increasing /x, and this tendency stops at about 
fi = 530MeV, which is related to the momentum cut-off A. 

Fig. |] is the phase diagram for color superconductivity when electric and color charge 
neutrality condition is considered, the quark mass m (squares), diquark condensate A (tri- 
angles) and the electron's chemical potential /i e (circles) are plotted as functions of quark's 
chemical potential /i. It is found that when chiral symmetry restores, the color super- 
conducting phase appears. In the color superconducting phase, it is found that both the 
diquark condensate A and the electron's chemical potential first increase with increasing /i, 
and reaches their maximum at about fi = 475 Me I 7 , then decrease rapidly with increasing 
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H, and the diquark condensate disappears at about /i = 535MeV. 

In Fig. [5], we show the chemical potential /igc as a function of /i, which ensures the 
system having zero net color charge T 8 . It is found that fi 8c can be negative and positive, 
but its value is very small, only about several MeV. It means that for the same flavor quark 
u or d, the difference of Fermi momenta between the third colored quark and the first two 
colored quarks can be neglected. 

In order to explicitly see what happens to a color superconducting phase when charge 
neutrality condition is added, we compare the electron's chemical potential in the neutral 
normal quark matter [i e (light circles) and in the neutral superconducting phase /i e (solid 
circles) in Fig. ||, and compare the diquark condensate A (light squares) and A (solid 
circles) in the electric charged and neutral superconducting phases in Fig. [7[ 

From Fig. ^|, we find that the electron has a larger Fermi surface in the neutral super- 
conducting phase than that in the neutral normal quark matter. This is because in the 
superconducting phase, the isospin is mainly carried by unpaired quarks only, and \x e has to 
be larger in order to get the same isospin density. \i e is equal to /i° at about /i = 535MeV, 
at which the diquark condensate disappears. We seperate /i e into two parts as fi e = fi^ + Sfig, 
where <5/i e reflects the effect induced by diquark condensate, comparing b\x e (light squares) 
and A (solid squares) in Fig. |6], we can see that <5/i e increases when A increases, and 
decreases when A decreases, and when A = 0, <5/i e = 0, too. 

From Fig. |7|, we can see that the diquark condensate A in the neutral superconducting 
phase has a much smaller value than A in the charged superconducting phase, it means 
that the difference of the Fermi surfaces of the two pairing quarks reduces the magnitude of 
the diquark condensate A. 

From Fig. [5] and Fig. |7|, we can see a distinguished characteristic in the neutral super- 
conducting phase, i.e., both the electon's chemical potential /i e and the magnitude of the 
diquark condensate A first increase with increasing then at about \i = 475MeV, decrease 
with increasing /i. Now we try to understand this phenomena. 

The magnitude of the diquark condensate is not only affected by /x, but also by 5[i = /x e /2, 
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which describes the difference of Fermi surface between the two paired quarks and reduces 
the diquark condensate as seen in Fig. [7]. The diquark condensate first increases with 
increasing fi, which is the result of the increasing density of states. At about fi = A75MeV, 
where fid = fi + l/3fi e is about 530MeV, the diquark condensate starts to be affected by the 
momentum cut-off A. Like that in the charged color superconducting phase, the diquark 
condensate stops increasing with fi, so does fi e . However, when fi > 475MeV, fi° e keeps 
increasing with fi, which reduces the diquark condensate largely. This is why we see A 
decreases with increasing fi. Finally at a certain chemical potential fi = 535MeV, where 
fid is about 5S0MeV, the diquark condenate disappears. Because Sfi e decreases with the 
decrease of A, the total fi e decreases in the chemical potential regime fi > 475MeV. 

At last, in Fig|| we show the thermodynamic potentials in the four different cases as 
functions of chemical potential fi, the light squares and circles are for the cases with and 
without charge neutrality and without diquark condensate, the dark circles and squares are 
for the cases with and without charge neutrality and with diquark condensate. It can be 
seen that in the chiral symmetric phase, the charged superconducting phase has the lowest 
thermodynamic potential, and the charged normal quark matter has the second lowest Q. 
The neutral superconducting phase has a little bit lower Q than that of the neutral quark 
matter in the chemical potential region fi < 535MeV, and at fi — 535MeV, the two fls 
coincide with each other. Therefore, in the chiral symmetric phase, the stable state of the 
neutral system is the color superconducting phase for < fi < 535MeV. At fi = 535MeV, 
the stable phase is the nomoral neutral quark matter. 

V. CONCLUSIONS 

We invetigated the effect of charge neutrality on a two flavor quark system. It has been 
found tha the BCS diquark pair has been largely reduced by the large Fermi momentum 
of electron, the diquark condensate first increases with the increase of quark's chemical 
potential fi, then decreases rapidly and diappears at about fi = 535MeV, at which the 



thermodynamic potential equals to that in the neutral normal quark matter. In the chem- 
ical potential region 330MeV < /i < 535MeV, the stable neutral system is in the color 
superconducting phase if no strange quark involved. 

As we mentioned in the introduction, we did not consider the strange quark in the system, 
which should be considered and has been investigated self-consistently in |J . Comparing our 
results with their results about the diquark gap A, the electron's chemical potential \i e and 
the chemical potential fi$ c in 2SC, it can be found that the main difference lies in the chemical 
potential region /i > AOOMeV. In [§], the diquark condensate in the neutral superconducting 
phase is largely reduced only in a small chemical potential region 370MeV < fi < AOOMeV. 
\i e first increases with \i then begins to decrease at about \i = AOOMeV. As for /i 8c , the 
tendency also becomes different from our results at about fi = AOOMeV. The reason lies in 
that, in ||, the strange quark involves in the system when [i > AOOMeV. This shows that 
the effect of the charge neutrality condition on three- and two- flavor quark system is quite 
different. 
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FIG. 1. The quark mass m as a function of chemical potential fi in the case of no color super- 
conducting phase and no charge neutrality condition. 
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FIG. 2. The quark mass m (squares) and the electron's chemical potential fi® (circles) as func- 
tions of chemical potential fi in the case of neutral normal quark matter. 
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FIG. 3. The quark mass m (squares) and the diquark condensate Ao (circles) as functions of 
chemical potential [i in the case of no charge neutrality on color superconducting phase. 
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FIG. 4. The quark mass m (squares), the diquark condensate A (triangles) and the electron's 
chemical potential [i e (circles) as functions of chemical potential fi in the case of considering charge 
neutrality on color superconducting phase. 
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FIG. 5. The figc as a function of chemical potential \x in the case of considering charge neutrality 
on color superconducting phase. 
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FIG. 6. The electron's chemical potential in the neutral normal quark matter ^ (light circles) 
and in the neutral superconducting phase fj, e (solid circles) as functions of chemical potential fi, 
and 5 ii e = ji e — (light squares) comparing with A (solid squares) as functions of jjl. 
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FIG. 7. The diquark condensate in the electric charged Ao (light squares) and in the neutral A 
(solid circles) superconducting phase as a function of chemical potential \x. 
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FIG. 8. The thermodynamic potentials as functions of chemical potential fi, the light squares 
and circles are for the cases with and without charge neutrality and without diquark condensate, 
the dark circles and squares are for the cases with and without charge neutrality and with diquark 
condensate. 
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